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Abstract 

We have studied the dissociation of quarkonia states in a deconhned medium of quarks and gluons 
at the large baryon chemical potential and small temperature region. The aim of this study is to 
probe the dense baryonic medium expected to be produced at the Facility for Anti-proton and Ion 
Research (FAIR), GSI Darmstadt. This is done by correcting both the short- and long-distance 
terms of the Cornell potential by a dielectric function, embodying the effects of deconhned quarks 
and gluons, at finite baryon chemical potential and temperature. It is found that J/ip is dissociated 
approximately at 1.4/x c in the temperature range 20-50 MeV, which can indirectly help to locate 
the point on the QCD phase diagram at the large chemical potential and low-temperature zone. 
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1 Introduction 

Relativistic heavy-ion collisions provide the unique opportunity of creating a hot and dense nuclear 
matter in the laboratory. According to statistical quantum chromodynamics (QCD), nuclear mat¬ 
ter may undergo a color deconhned partonic phase, the quark gluon plasma (QGP), at sufficiently 
high temperature and/or density. Over the past few decades, strenuous efforts have been made 
to devise clean and experimentally viable signals that can unambiguously identify the existence 
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of QCD phase transition and trace out its signatures. Charmonium (a bound state of charm and 
anti-charm quarks) suppression had been predicted as a signature for the deconfinement transition 
for a long time [1], where the Debye screening of the colored partons was the dominant mechanism 
of suppression. Experimental investigations revealed a noticeable reduction of the charmonium 
production in proton-nucleus (p-A) collisions compared to the hadronic collisions (scaled with 
binary nucleon-nucleon (p-p) collisions) [2, 3, 4, 5]. A thorough understanding of this normal 
nuclear suppression is essential in order to establish a robust baseline reference, with respect to 
the anomalous suppression [6, 7], pertaining to the formation of a deconfined partonic medium 
observed at Relativistic Heavy Ion Collision (RHIC) and Large Hadron Collider (LHC) energies 
that mainly scan the high temperature and almost vanishing baryon density region of the QCD 
phase diagram [8, 9]. Depending upon the magnitude of energy deposited in the center-of-mass 
frame of colliding nuclei, there are in general two scenarios in heavy-ion collision. The first one 
is possibility of the complete stopping of the colliding nuclei if energy density e ~ 5 GeV/fm 3 , 
creating sufficiently high baryon density (parton density) to produce QGP. Full transparency of 
colliding nuclei, the other scenario, takes place if e ~ 100 GeV/fm 3 creating a favorable envi¬ 
ronment to produce QGP by vacuum polarization and leading to vanishing chemical potential. 
Exploration of high baryon densities and the moderate-temperature region of the QCD phase dia¬ 
gram is possible with the upcoming compressed baryonic matter (CBM) experiment at the Facility 
for Anti-proton and Ion Research (FAIR). At FAIR, the light and heavy ions can be collided in 
the beam energy range E ^ = 10-40 AGeV to create an environment favorable to investigation 
of extremely dense nuclear matter in the laboratory through the measurements of bulk and rare 
probes. Model calculations based on transport and hydro dynamical equations [10] predict that 
the highest net baryon densities ( ps ) produced in the center of collision is ~ 6 to 12 times the 
density of normal nuclear matter for the most central collision (b=0). 

The FAIR will open avenues to investigate some of the fundamental yet enigmatic issues of 
strong interaction thermodynamics at large baryon density. Some of the key issues in this region 
are the study of hadronic properties in dense nuclear matter, the deconfinement phase transition 
from hadronic to quark-gluon matter driven by high baryon densities, and the nuclear equation of 
state (EOS) at high baryon densities. Thus the behavior of QCD at high baryon density and mod¬ 
erate to low temperatures is interesting and has potential applications to cosmology, astrophysics 
of neutron stars, and heavy-ion collisions. Although the motive of FAIR is to produce baryon-rich 
QGP, it is uncertain that at what point on the phase diagram along the baryon density ( Pb) axis 
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QGP can be produced. 


Since the relative velocity of heavy-quark-bound sate (vqq 1) is very small, in-medium 

dynamics of QQ bound states has been extensively studied with phenomenological potential mod¬ 
els [11], where the temperature-dependent potential carries all medium effects with nonperturba- 
tive terms borrowed from lattice simulations. Although the derivation of such models from QCD 
is not established, free energies and other quantities [12, 13] from lattice calculations obtained 
using the correlation functions of the Polyakov loop are often taken as input for the potential. 
These quantities have been thought to be related to the color-singlet and color-octet heavy quark 
potentials at finite temperature [13, 14] however, a precise answer is still missing [15]. It was 
until recently that the smallness of the velocity ( Vqq ~ a s < 1) has opened up the possibility 
of studying heavy quark bound sate at finite temperature, where both the non-relativistic scales, 
viz. the heavy quark mass (mg), the momentum exchange (mQfg), the binding energy (mgWg) 
etc. and the thermal scales, viz. T , gT , g 2 T are are hierarchically ordered: By exploiting the 
aforesaid separation of scales, a sequence of low-energy effective held theories (EFTs) [16, 17] 
viz. non-relativistic QCD (NRQCD), potential non-relativistic QCD (pNRQCD) etc. have been 
synthesized by integrating out the successive energy scales. In this context, pNRQCD and its 
thermal version are of particular interest because it describes the quarkonium dynamics through 
potentials and low-energy interactions [18, 19]. The color-screening phenomenon [18] and thermal 
width are thought to be an outcome of thermal corrections to the real and imaginary parts of the 
color-singlet potential, respectively. Out of various mechanisms that might be contributing to the 
imaginary part of the potential, Landau-damping [20] is the dominant one, which is responsible 
for quarkonium dissociation at weak coupling. It makes the quarkonium to be dissociated even at 
temperatures, where the possibility of color screening is negligible. Moreover, lattice studies have 
also shown that even at strong coupling [21, 22, 23] the potential may have a sizeable imaginary 
part and such contributions are related to quarkonium decay processes in the plasma. 


The works referred above are restricted to finite temperature only. However, the color screening 
effect at finite temperature and finite density has been studied in thermo-field dynamics approach 
to calculate the Debye screening mass, mfl(T, /i), where authors in Refs. [24] and [25] have used the 
phenomenological potential model [11] and error function-type confined force with color screened 
Coulomb-type potential, respectively. Applying the effective perturbation theory for gauge theo- 
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ries at finite temperature [26], authors in Ref. [27] extended it to the finite chemical potential for 
studying the collisional energy-loss of heavy quark in QGP. Photon production at finite chemical 
potential in QCD plasma at leading-order in strong coupling has been computed in [28]. Further¬ 
more the dissipative hydrodynamic effects on QGP at finite density has been studied by developing 
a causal dissipative hydrodynamic model at finite baryon density for RHIC and LHC energies to 
study the net-baryon rapidity distribution [29]. Although, finite density lattice QCD calculations 
are seriously affected by the sign problem, recently color screening in heavy quark potential at 
finite density with Wilson fermions has been studied in lattice QCD in both real and imaginary 
chemical potential regions [30]. 


In the present work we have studied quarkonium dissociation in the quark matter of high 
baryon density and low temperature. For that we have first obtained the heavy quark potential 
by correcting both perturbative and nonperturbative terms in the Cornell potential through the 
dielectric function at finite temperature and chemical potential and then calculate the binding 
energy of the quarkonium states. The rest of the work is organized as follows: We first discuss the 
heavy quark potential at finite chemical potential in Sec.II. Then we focus on the estimation of 
the binding energy of quarkonium at finite chemical potential and thereby obtain the dissociation 
chemical potential for quarkonium ground states in Sec.III. Finally we conclude in Sec.IV. 

2 Heavy quark potential at finite chemical potential and 
running coupling constant 

The interaction potential between Q and Q is subjected to modification when the pair (QQ) is 
placed in a hot and dense QCD medium, and it plays an important role in understanding the 
status of QQ bound states in such medium. This issue has been well taken up by various authors 
and is reported in several reviews [31, 32], by dealing on both fronts, the phenomenology [11] as 
well as the lattice QCD [21, 22, 23, 33]. All these studies assume that the phase transition from a 
hadronic matter to a QGP phase melts the string, and thereby the string tension vanishes at the 
transition point. However, lattice results [34] indicate that there is no genuine phase transition 
at vanishing baryon density, it is rather a crossover, and there may not be immediate melting of 
string at the deconfinement temperature. Thus one should consider its effect on the potential even 
above the transition temperature, in addition to the coulomb term. 
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The large mass of heavy quark meets the requirements: (i) mg A qcd and (ii) T -C tuq 
for the description of the interactions between QQ at finite temperature and density, ft is thus 
possible to obtain the quantum mechanical potential V(r,T), by correcting short as well as long¬ 
distance parts of the QQ potential through dielectric permittivity, e(k), embodying the effects of 
deconhned medium: 


V{r,T,ii) 



( 1 ) 


The term V(k) in above equation is the Fourier Transform (FT) of the potential. However, 
obtaining the FT of the linear term is a tricky part and needs regulation, which is carried out by 
multiplying an exponential damping factor and is then switched off after the evaluation of FT. 
We have regulated both terms by the same screening scale 3 and the FT of the Cornell potential 
is thus obtained as [35] 


V(fc) 


Vvl 

k 2 



( 2 ) 


where a and a are the coupling and string tension, respectively. We have taken a as function of 
chemical potential (fi q ) and temperature up to two loops [36] as 
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where is the QCD scale-fixing parameter which characterizes the strength of the interaction. To 
obtain the dielectric permittivity, e(k), one has to evaluate the self-energies and the corresponding 
static propagators in weak coupling Hard Thermal Loop (HTL) approximation. In the real¬ 
time formalism using Keldysh representation, the retarded (R), advanced (A) and symmetric 
(F) propagators are written as the linear combination of the components of the (2 x 2) matrix 
propagator: 


D° r = D° u 


D° 

-^12 


nO _ pO 


D° 

^21 


n° — no 

u F — u ll 


D° 
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( 4 ) 


where only the symmetric component involves the distribution functions and is of particular advan¬ 
tage for the HTL diagrams, where the terms containing distribution functions dominate. Similar 
relations hold good for the retarded (n«), advanced (n^) and symmetric (n^) self energies. The 

3 However, in the framework of classical Debye-Hiickel theory, different screening functions were employed, viz. 
f c and f s for the Coulomb and string terms, respectively to obtain the free energy [37] whereas in Ref. [38], different 
scales for the Coulomb and linear pieces have been employed through a dimension-two gluon condensate. 
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resummation of the propagators can be obtained using the Dyson-Schwinger equation 


(5) 

( 6 ) 


-Dr,a — D° ra + D° R A U RtA D RtA , 

D f = D° f + D° r U r D f + D° f TI a D a + D° r U f D a . 


For the static potential, we need only the temporal component (“00” = L) of the propagator, 
whose evaluation is easier in the Coulomb gauge. Thus the above resummation (5) can be recast 
through its temporal component as 


Dr,a = D 


L(0) 

R,A 


d%°) n L 


R,A ll R,A j ^R,A 


D\ 


(7) 


The leading contribution to temporal component of the retarded, advanced and symmetric gluon 
self-energy in the HTL-approximation can be written as 



m 2 D 


hi„ 

2k 


k 0 + k ±ie 
k 0 — k ±ie 



( 8 ) 


and 


H F {k) = -27rim 2 D ^<d(k 2 - k 0 2 ) , (9) 

respectively. Thus the gluon self-energy is composed of real and imaginary parts which are respon¬ 
sible for the Debye screening and the Landau damping, respectively. The real part of potential 
can be obtained from the retarded (or advanced) self energy with the prescriptions +ie (— ie), 
respectively. 

The term m D in above equation is known as electric screening mass (Debye screening mass) which 
encodes the medium effects in terms of temperature and baryon chemical potential. The properties 
of QCD at finite temperature and density are studied in [39] and the electric screening mass can 
be obtained in the static infra-red limit of “00” component of gluon self energy II 00 (g 0 = 0, q—> 0). 
Taking the dynamical quark mass m q = 0, the electric screening mass at finite temperature and van¬ 
ishing baryon chemical potential {/i = /if, = 0) based on perturbation theory in high temperature 
limit is [40], 

n 00 (<?o = 0,q 0) = ^ ^ gT, (10) 

where N c and Nf are number of colors and flavours of color group SU (N c ) respectively. For a hot 
and dense plasma, the resummation technique to the perturbation theory proposed by [26] works 
well in momentum scale gT. The leading-order Debye screening mass from HTL [27, 41, 42, 43] 
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at finite temperature and chemical potential is given by 

m 2 D (T,fi q ) =g 2 {T)T 2 j ’ ( n ) 

where the quark chemical potential is related to the baryon chemical potential by /j, q = //&/3. 
Later the screening mass at finite temperature and baryon density has also been studied on the 
lattice through the Taylor expansion method [44], which coincides with the above equation. For 
zero chemical potential (/i q = 0), the leading-order Debye mass (11) reduces to 

m D (T)=g(T)T ] /^ + ^-. (12) 


The next-to-leading contribution to Debye mass (5mi)(T, fi q )) comes from the resumed one- 
loop gluon diagrams [26]. However, we restrict ourselves only up to the leading-order term in the 
Debye mass. The variations of the Debye mass in leading-order with the chemical potential and 
temperature are shown in Fig. 1. The left panel of Fig.1 shows that the Debye mass increases with 
the chemical potential and is more for higher temperature. This implies that the contribution to 
the Debye mass at low temperature and high baryon chemical potential mainly comes from the 
chemical potential, which can be seen in the right panel of Fig.l. For low temperature region (20- 
50 MeV) the increase in Debye mass is negligible but increases with / 15 . The temporal component 
of the retarded (or advanced) propagator in the static limit [45] is 


D° r ° a ( 0,k) 


1 

(P + m^T,/!))’ 


(13) 


Hence we can now obtain the dielectric permittivity from the static limit of the “00”-component 
of gluon propagator 


e~\k) 


— lim k 2 D^(uj, k) 
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Thus the medium modified potential is obtained by substituting Eqs. (2) and (14) into Eq.(l) 


V(t ,T,/i) = 


d 3 k 
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Solving the integral in above equation, one can to obtain the medium modified potential at finite 
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Figure 1: Variation of the screening mass (mu) against the baryon chemical potential (fib) at 
different values of temperature (left panel). Same against the temperature at different values of 
baryon chemical potential (right panel). 


density and temperature (with r = rmr,(T, /i)) as: 

V(r, T, fi) = - . 1 + l\ - am D (T,n) (%r + l) , (16) 

V r J V r J 

where the chemical potential and temperature dependencies are introduced through the Debye 
mass, mo(T, fi). The potential has a long range Coulombic tail in addition to the standard 

Yukawa term and the constant terms are introduced to yield the correct limit of V (r, T, /i) as 

T —* 0, /i —* 0 (it reduces to the Cornell potential). 


Recently it is known that there exists imaginary part of the potential [16, 17, 21, 46] which 
originates from the static limit of symmetric self energy (9) and plays an important role in weak¬ 
ening the bound state peak that leads to a finite width (T) for the resonance peak in the spectral 
function. The imaginary part of the temporal component of symmetric propagator gives the 
corresponding dielectric function [45] 


e 1 (0, fc) 


— TiTm 2 D (T, fi)k 2 
k(k 2 + m 2 D (T, /i)) 2 ’ 


(17) 


which, in turn, gives the imaginary part of QQ potential: 
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where the first (second) term comes from contribution due to the Coulomb (linear) term in the 
leading-order [47, 48, 49] with the following functions, 0 O if) and ipo (f), respectively: 


Mr) 

Mr) 


-■^-(-4 + 37 E + 31ogf) 

—107 + 607 ^ + 60 log(r) 5 

6 3600 1 J 


(19) 

( 20 ) 


Thus in leading-logarithmic order, the imaginary component of potential becomes 


QV(r,T, M ) = -T 


ar 


or 4 \ 1 

+ 30m 2 D (T, n)) ° s( F' 


( 21 ) 


which shows that imaginary part vanishes for small distances [45]. However, its magnitude is larger 
than the case where only the Coulombic term is considered [47, 50, 51], due to the additional linear 
(string) term. 

Now with the imaginary potential (21), the width for a particular resonance state can be calculated 
as: 


rp »/2 



Tat 4 


30 rrijy(T, /i) 

arriQ 


log(l/r) 


m D log 


2 m D (T, n)' 


( 22 ) 


Here for simplicity, we have taken T(r) as the Coulombic wave function, similar to the ground 
state wave function in hydrogen atom. This shows that T is function of both chemical potential 
and temperature through the screening mass (Debye mass). However, the dependences on tem¬ 
perature and chemical potential are opposite, namely V increases with temperature but decreases 
with chemical potential. That is why T does not have much role in dissociation of quarkonia states 
unlike at high temperature. This will be discussed in the next section. 

It is worth to note that the potential in a hot QCD medium is not the same as the lattice 
parametrized heavy quark free-energy in the deconhned phase which is basically a screened 
Coulomb [52] because one-dimensional Fourier transform of the Cornell potential in the medium 
yields the similar form as used in the lattice QCD to study the quarkonium properties which 
assumes the one-dimensional color flux tube structure [53]. However, at finite temperature and 
density that may not be the case since the flux tube structure may expand in more dimensions [52], 
hence the three-dimensional form of the medium modified Cornell potential may therefore be the 
better option. The potential (16) leads to an analytically solvable Coulomb potential if one neglect 
Yukawa term in the limit r » 1 /me, and the product amp will be much greater than 2a /mjj 
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Figure 2: Variation of medium modified QQ potential in a hot and dense QCD medium at given 
low temperature with different baryon chemical potential /x;, (left panel). Same for different tem¬ 
peratures at given fib (right panel). 


for large values of chemical potentials. 

V ( r , T,fi )-- am D (T, fi) (23) 

To see the effects of a hot and dense medium on the QQ potential, we have evaluated the potential 
at large chemical potential, viz. at 600 MeV, 800 MeV and 1000 MeV in the low temperature 
range 20-60 MeV in Fig.2. We have found that the potential rises with the interquark distance, 
r, but slower than linearly, i.e. long-range QCD force becomes short-range and at some large 
enough r, it simply flattens out due to the breaking of QCD string, i.e. a heavy meson splits into 
two heavy-light mesons. This observation agrees with the lattice calculations [54, 55]. Thus, the 
deconfinement is reflected in the large-distance behaviour of heavy quark potential even at large 
chemical potential. That is why the in-medium behaviour of heavy quark bound states is used to 
probe the state of matter in QCD thermodynamics at finite density and/or temperature. 

3 Binding energy and dissociation chemical potential 

The medium-modified effective potential obtained in the previous section has been employed to 
study the in-medium properties of the heavy quark bound states such as, binding energies and 
dissociation chemical potentials. We studied the binding energy and the dissociation chemical 
potential for the ground states of cc. In order to understand bound state properties of quarkonium 
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states in hot and dense QCD medium, one need to solve the Schrodinger equation using the QQ 
potential. It can be seen from Eq.(16) that in the short-distance limit, the vacuum contribution 
of QQ potential supersedes over the in-medium contribution whereas the in-medium contribution 
affects the potential in the long-distance limit. As we have already seen the potential (23) in 
long-distance ( rm,D 1) limit and for large values of chemical potential (or temperature) is 
Coulomb-like potential after identifying 2cr/m 2 D with the square of color charge, g 2 , so the energy 
eigenvalues are read as in hydrogen atom problem: 


P *»i 

& bin — 


m 


2iyv +Qm ‘ ,(3>) 


n = 1,2 ■■■ 


(24) 


However, in the intermediate-distance (rmi)(T, fi) ~ 1) scale, the interaction becomes non¬ 
trivial and the potential does not look simpler in contrast to the asymptotic limits, thus the full 
potential (16) in general needs to be dealt numerically to obtain the binding energies. The matrix 
method serves the purpose and the stationary Schrodinger equation can be solved in a matrix 
form through a discrete basis, instead of the continuous real-space position basis spanned by the 
states |T’). Subdividing the potential V (16) into N discrete wells having potentials V\ through 
Vjv+ 2 , such that V = Vi for Xj_i < x < Xi\ i = 2, 3,..., (N + 1) is the i th boundary potential. For 
existence for a bound state, there must be an exponentially decaying wave function in the region 
x > x^v+i as x —> oo and the wave function in this region reads: 


Tjv+sOr) = P E exp[—-7 iV+2 (x - xjv+i)] + Q E exp[7 JV+2 (x - xjv+i)], ( 25) 

where, P E = \{A N+2 - B N+2 ), Q E = \{A N+2 + B N+2 ) and, -f N+2 = \/2fi r (V N+2 - E). The energy 
eigenvalues can then be identified with the zero’s of Qe- 

Thus, matrix method gives the binding energies as a function of chemical potential in a small 
temperature range (20 - 40 MeV), which are shown in Fig. 3. It is found that the binding energy 
decreases with the increase of chemical potential. This can be understood by the fact that the 
screening gets stronger with the increase of the baryon chemical potential and hence the inter¬ 
quark potential becomes weaker compared to fi = 0. Thus the study of the binding energies at 
finite temperature and chemical potential provides a information about the dissociation pattern 
of quarkonium states in the region of high baryon chemical potential and low temperature. 

Now one can obtain the dissociation chemical potential of quarkonia when their binding energies 
in the dense baryonic medium are of the order of the baryon chemical potential. We found that 
with the rise in temperature, the dissociation chemical potential deceases (Table 1) slightly, for 
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Figure 3: Variation of binding energy of J/rj) state in a hot/dense QCD medium with baryon 
chemical potential at different values of temperature. 


T (MeV) 

Pd (MeV) 

/ ftc 

20 

1683 

1.42 

30 

1659 

1.40 

40 

1611 

1.38 

50 

1560 

1.35 


Table 1: Dissociation chemical potential for J/rj) at different temperatures. 


example, at temperature T = 20 MeV, the ■J/'tp will be dissolved at 1.42 p c whereas for T=30 Mev, 
Pd becomes 1.40 /i c 4 etc. This observation can be used to locate the point on QCD phase diagram 
in the high-baryon density limit where the confined nuclear matter cross to the phase boundary 
and becomes baryon-rich deconhned quark matter Recently the physics of dissociation has been 
further complicated by arguing that the potential should have an imaginary part as well, which, 
like a width, facilitates the dissociation of quarkonia further. Thus, nowadays the dissociation 
temperatures are obtained from the width, T, where the dissociation is expected to occur when 
binding energy decreases with the chemical potential and becomes ~ T/2 [56, 57]. We have found 
that the dissociation chemical potential estimated from the above criterion is almost same to the 
values obtained from the binding energy criterion (Table 1). 

4 The values of T c and p, c used in this work are taken from the Ref. [58], viz. the critical baryon chemical 
potentials (p c ) are 1189 MeV, 1183 MeV, 1169 MeV and 1154 MeV for the critical temperatures (T c ) 20 MeV, 30 
MeV, 40 MeV and 50 MeV, respectively. 
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4 Conclusion 


We have studied the quarkonium dissociation in hot and dense QCD medium by correcting both 
the perturbative and nonperturbative terms of the Cornell potential through the dielectric per¬ 
mittivity Thereafter the chemical potential- and temperature-dependent potential is plugged into 
Schrodinger equation to study the properties of quarkonia states in high baryon density and low 
temperature region. 

ft is noticed that the screening mass increases dominantly with the chemical potential rather 
than with the temperature in the high baryon density and low temperature region. As a con¬ 
sequence the binding energies decrease with the chemical potential. The dissociation chemical 
potentials of J/ijj at T= 40 and 50 MeV is found to be 1611 MeV and 1560 MeV, respectively. 
Finally, we conclude that the dissociation of quarkonia states in large baryon density may help to 
trace out the point on the QCD phase diagram at which the baryon rich QGP (which is expected 
at FAIR energies) will be created. 
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